Interstellar Dust Emission as a CMBR Foreground 
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Introduction 



Abstract. 

I This chapter discusses Galactic dust and how its thermal emission 

confuses CMBR measurements. Interstellar dust grains are composed of 
many differing species, and observational evidence has only begun to dis- 
entagle their composition and sizes. Fortunately, their far-IR/millimeter 
I emission is less complex. We describe how a very simple two-species 

' model can describe the emission at 200 < < 3000 GHz to high preci- 

I sion. At lower frequencies, other non-thermal processes may dominate 

' the emission from dust. 

o 

' The different sources of radiation in the night sky from the near-infrared through 

, radio are shown in Figure |l[ These are to be compared with the CMBR at a 

temperature of 2.73 K, and the fluctuations in the CMBR which have an am- 
plitude lower by ~ 10"^ (dashed line). On the Wien tail of the CMBR, thermal 
^ ' emission from Galactic dust is the main source of confusion. This emission is 

now well-understood, thanks to the COBE/FIRAS experiment (Fixsen et al. 
1996). We can predict the thermal emission at all frequencies ^ 3000 GHz 
to a precision of ~ 10% (Finkbeiner, Schlegel & Davis 1999). This work builds 
upon the maps of dust column density produced by Schlegel, Finkbeiner, &; 
Davis (1998; hereafter SFD98). 

At frequencies below ~ 100 GHz, there are more possible sources of con- 
fusion with the CMBR. Galactic dust may be one such source if it is rotating 
supra-thermally, as has been proposed by Draine & Lazarian (1998). Other pos- 
sible sources are synchrotron emission, free-free emission from ionized hydrogen, 
and radio emission from extragalactic sources. We briefly discuss all as possible 
sources for anomalous emission seen at ^ 100 GHz by CMBR experiments. 

1.1. Observational Evidence 

Dark clouds are apparent to any observer who looks at the Milky Way from 
a dark site. These dark regions are especially conspicuous in the Southern 
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Figure 1. Contributions to the infrared/millimeter /radio sky above 
the Earth's atmosphere 

hemisphere where one sees the utter darkness of the Coal Sack and dark bands 
across Centaurus and Sagittarius. Wilham Herschel regarded these regions as 
curious voids in the distribution of stars. In the early 20th century, detailed 
studies by Trumpler convincingly demonstrated that interstellar dust was re- 
sponsible. Starlight passing through a dusty medium is subjected to both ab- 
sorption and scattering. The combination of these effects is termed extinction, 
denoted A{V) for the magnitudes of visual-band extinction. Both processes are 
wavelength-dependent, resulting in an apparent reddening of the light. Red- 
dening is usually parameterized by its differential extinction between B- and 
y-band, E{B - V) = A{B) - A{V). 

Advances in infrared astronomy since 1980 have dramatically improved our 
understanding of interstellar dust. Dust scatters some light in the optical, as 
can be seen in the Pleiades star cluster. In our solar system, a thick disk of 
dust in the zodiacal plane scatters some sunlight. However, more energy is 
absorbed and then re-radiated in the infrared. The zodiacal dust re-radiates 
most of its energy at A ;^ lO^m. The energetics of this zodiacal dust would 
be uninteresting except that it dominates the sky (above the atmosphere) at 
wavelengths ~ 25/im. This zodiacal emission has been extensively modeled by 
Kelsall et al. 1998. Interstellar dust is heated to a lower temperature, radiating 
its energy at A ^ lOO^m. Approximately 20% of the energy from starlight in 
our Galaxy is re-processed this way by interstellar dust. 
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Dust Composition The diffuse ISM is known to contain many different types 
of molecules and dust grains whose optical properties are dissimilar. Identifica- 
tions for some of these constituents have been made from UV absorption features 
and mid-IR emission features. The ubiquitous 2175A bump (in absorption from 
background stars) is most widely attributed to small (< lOOA) graphite parti- 
cles (c./., Draine 1989). Carbon in another form, that of polycyclic aromatic 
hydrocarbons (PAHs), is probably responsible for strong emission features at 
3.3, 6.2, 7.7, 8.6 and 11.3;um (e.g., Allamandola, Ticlcns, & Barker 1985; Desert 
et al. 1990). A review by Henning & Salama (1998) contains the latest infor- 
mation on carbon-based grains. Other mid-IR emission features are attributed 
to silicates and troilites. For a coherent and only somewhat dated review, see 
Whittet (1992). 

Grain Sizes Along with the composition of dust grains, their size distribution 
also impacts their extinction and emission properties. The UV/optical extinc- 
tion, ^(A), is neatly described by a one-parameter family of curves (Cardelli, 
Clayton, & Mathis 1989). That one parameter is usually chosen as an observable 
ratio, 

_ AjV) 

- A{V)-A{By 

Clayton & Cardelli (1988) have measured the extinction curve for 76 lines-of- 
sight in our Galaxy. Most diffuse clouds favor an Ry = 3.1 extinction curve. 
However, there are variations in Ry from 2.6 to 5.5, where the larger values 
correspond to some dense cloud environments. Changes in grain size distribu- 
tion is presumed responsible for these variations. On average, interstellar dust 
grains are thought to follow a power law distribution of grain sizes dn/da oc a~^-^ 
(Mathis, Rumpl, & Nordsieck 1977) from about 0.0005/im (the size at which ab- 
sorption of a single photon can sublimate mass away from the grain) to 0.25/xm. 
The large-size cutoff at 0.25/Ltm is necessary to be consistent with the measure- 
ments of Ry (Kim, Martin & Hendry 1994). Larger grains would be relatively 
gray absorbers, increasing the value of Ry beyond that observed. So Ry pro- 
vides crucial information on the grain size distribution, and may be found in the 
future to correlate with some aspect of the far-IR emission from dust. 

Polarization Effects Besides attenuating starlight, dust grains are known to 
polarize that light. To date, polarization of optical light has been measured 
towards ~ 10,000 stars (Mathewson & Ford 1970; Korhonen & Reiz 1986), 
some high above the Galactic plane (Berdyugin & Teerikorpi 1997). Regions of 
higher polarization are correlated with regions of higher extinction (e.g., Larson, 
Whittet, &; Hough 1996). One expects polarization only if the dust grains are 
non-spherical and have a preferential axis of alignment. For this reason, polar- 
ization studies have been used to argue for coherent, large-scale magnetic fields 
in the Galaxy {e.g., Heiles 1996). The magnetic moments of the dust grains will 
preferentially align with the ambient magnetic field. As these aligned grains ro- 
tate, they emit polarized radiation that may interfere with attempts to measure 
the polarization of the CMBR anisotropy. 

Dust Distribution The morphology of the interstellar dust is quite complex, 
consisting of dark clouds and long filamentary structures arching across the sky. 
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The diffuse, wispy parts of this emission are known as cirrus, to differentiate 
them from more compact features (often associated with individual objects). 
Not surprisingly, the dust emission is concentrated in the Galactic plane, but 
significant cirrus is observed in every part of the sky at lOO/um (Low et al. 
1984). SFD98 have computed the power spectrum of this dust at latitudes 
\b\ > 45°, finding a reasonable fit to P{k) oc k~^-^. However, as this cirrus 
is in coherent structures, a 1-D power spectrum is not statistically sufficient 
to describe its distribution. This morphology is another obstacle to CMBR 
measurements since it is not at all reminiscent of gaussian random-fields that 
CMBR pundits usually model. However, this phase coherence provides avenues 
that may prove helpful in determining the quality of a dust foreground model 
as applied to real microwave data. For example, we can easily measure the 3- 
point function of the Galactic dust. If the CMBR is presumed to be gaussian, a 
non-zero 3-point function in a foreground-subtracted CMBR map will measure 
the level of residual contamination (Spergel 1998). 

1.2. Far-IR Emission 

In spite of the expected melange of dust grains, it was originally expected that in 
the far-IR/submillimeter bands, all dust would have similar optical properties. 
For example, Draine & Lee (1984) predicted v"^ emissivity for both silicate and 
graphite grains whereby the emission is described by 

7, = B,{Ty. (2) 

The theoretical motivation for this emissivity is based on the fact that the re- 
sponse of a damped harmonic oscillator to a driving force varies as in the 
low frequency limit (far below resonance). Even though microwave frequencies 
are below all of the phonon modes in the smallest crystals, this is not true for 
grains bigger than ~ 0.1/xm, and these bigger grains are expected to dominate 
the emission. More recent laboratory measurements (sec Pollack et al. 1994 
for a review) suggest that this is an oversimplication, with dift'crcnt species of 
grains having differing emissivity laws in the far-IR. In particular, glassy mate- 
rials (e.g. amorphous silicates) in some cases radiate more efficiently than their 
crystalline counterparts, and are more closely characterized by a i^^'^ emissivity 
(Agladze et al. 1996). As we show below, these silicates probably dominate the 
millimeter /microwave emission relevant to CMBR studies. 

Fortunately, complete knowledge of the size and composition of interstel- 
lar dust grains is not necessary to predict their impact upon CMBR studies. 
Since the long wavelength emissivity of a grain scales as its size a times the 
surface area, or volume a^, the larger grains dominate the submillimeter emis- 
sion (Draine Sz Anderson 1985) The very small grains (VSGs) are transiently 
heated to higher temperatures and dominate the dust emission at A ^ 100/um, 
but do not contribute significantly to submillimeter emission. The opacities of 
dust grains are exceedingly complicated in the UV/optical where they absorb 
energy, and slight changes in composition can significantly alter them. However, 
their far-IR opacities are not nearly as sensitive. Both laboratory and astronom- 
ical measurements are well-fit at A ;^ 1000//m {u ^ 300 GHz) using power-law 
emissivities with indices in the range from 1 to 3 (see Pollack et al. 1994). 
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2. Microwave Foregrounds 

Interstellar dust is only one of the many foregrounds that complicate measure- 
ment of the CMBR anisotropy. An analysis of the relative importance of the 
various foregrounds is given by Tegmark &; Efstathiou (1996). In this section, 
we provide a brief overveiw of the most important foregrounds. 

At low frequencies, the dominant diffuse Galactic emission is due to syn- 
chrotron radiation produced by relativistic electrons interacting with the Galac- 
tic magnetic field. A map at 408 MHz by Haslam et al. (1981) is the only 
full-sky synchrotron template in existence. Reich &; Reich (1986) cover the 
northern sky in the 1420 MHz continuum, and are complemented by the new 
Rhodes/HartRAO survey at 2326 MHz in the south (Jonas, Baart, & Nicolson 
1998). Where these maps all overlap, they allow examination of the synchrotron 
radiation at 3 frequencies ~ and indicate significant spectral index variation 
across the sky. The synchrotron radiation is usually characterized by a spec- 
tral index /? = —2.7 in units of brightness temperature, oi (3 = —0.7 in units of 
flux {e.g., MJy/sr). The typical high-latitude brightness temperature at 20 GHz 
predicted from these data is approximately 0.5 — 1 mK, rising to no more than 
20 mK in the Galactic plane. Extrapolating the power law to 100 GHz, the 
synchrotron emission would contribute at most 10 fiK to the emission at high 
latitude. However, the index (3 is seen to steepen to —2.9 or less even at low 
frequencies, due to a cutoff in the electron energy distribution (Jonas 1998). 
Therefore, synchrotron emission is unlikely to be an important emission compo- 
nent at high latitude at ^ 30 GHz. 

Another diffuse component that is not well constrained is free-free emis- 
sion from ionized gas. Except in Hll regions, free-free is lost in the signal of 
synchrotron emission at low frequencies {e.g., at 408 MHz). Although there is 
known to be ionized gas at high latitude from Ha emission (Haffner, Reynolds, 
& Tufte 1998), so far no convincing correlation has been seen between Ha data 
and ~ 20 GHz data. As more data of both types become available, we expect to 
see a correlation. Another possible emitter of free-free is the Local Hot Bubble 
(LHB) of ionized gas in which the Sun and many other stars reside. A map 
of this bubble was produced using the absorption of soft X-rays observed by 
ROSAT (Snowden et al. 1998). Because the gas in the LHB is very hot, the 
anticipated Ha emission per microwave emission is very low, so this will not 
appear in the Ha template. 

In addition to these diffuse components, a wealth of point sources are ap- 
parent at radio and far-IR frequencies. Several hundred thousand radio sources 
have been observed in the PMN (Griffith k Wright 1993) and NVSS (Condon 
et al. 1998) surveys. At 60^m, SFD98 identified ~ 20, 000 extragalactic sources 
and unresolved Galactic sources. These sources will contribute a great deal of 
non-gaussian contamination to the CMBR. The brightest ones are very easy to 
identify and remove. Preliminary work has shown that removing the brightest 
sources will be sufficient for power-spectrum analyses with the MAP experi- 
ment (Tegmark & de Oliveira-Costa 1998). Point source decrements due to the 
Sunyaev-Zeldovich effect will have also been seen, and their impact on CMBR 
observations has been noted (Aghanim et al. 1997). 

As we have emphasized, the dust emission described in the next section 
is not the only diffuse foreground, but it is the brightest at high frequencies 
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(z^ ^ 100 GHz). However, it is important for another reason. The excess dust- 
correlated emission observed in the DMR 31.5 and 53 GHz channels (Kogut et al. 
1996) and later by other experiments (e.g. Leitch et al. 1997; de Oliveira-Costa 
et al. 1998) is still unexplained. It may be caused by dust-correlated free-free 
emission. On the other hand, it may be the result of rapidly rotating dust grains 
(Draine & Lazarian 1998). Whatever this component is, it appears to correlate 
fairly well with dust, making the determination of a good dust template more 
important than thermal (vibrational) dust emission alone would indicate. 



3. Single Component Dust Model 

A simple but naive prediction for millimeter /microwave dust emission can be 
made from our previous work. SFD98 extensively studied the emission from 
dust in the regime lOO^um < A < 240//m. Assuming a emissivity model, the 
temperature of the dust was mapped with a resolution of 1.3° from the DIRBE 
lOO/im/240/im emission ratio. The lOO^m emission of the dust was mapped 
with a resolution of 6.1' by utilizing small-scale information from the IRAS Sky 
Survey Atlas (Wheelock et al. 1994). Emission at lower (millimeter /microwave) 
frequencies can be predicted by extrapolating the 100/um flux using this temper- 
ature fit. For each line-of-sight in the maps, the emission at frequency v can be 
expressed as 

I. = i^roJ,(a,r)/ioo5^|^ (3) 

where By{T) is the Planck function at temperature T, Iioo is the DIRBE- 
calibrated 100/xm map, KiQQ{a,T) is the color correction factor for the DIRBE 
100;um filter when observing a v'^ emissivity profile, and vq = 3000 GHz is the 
reference frequency corresponding to 100/im. 

The extrapolated, millimeter emission from dust as predicted from SFD98 
can be compared directly to the FIRAS measurements. We synthesize a broad- 
band FIRAS 500 GHz channel (to increase S/N) and show that the correlation 
with predictions from equation |3| is very good (Figure ^). For comparison, we 
plot the correlation with (a) Hi column density, and (b) DIRBE lOO/im flux. 
Both of these show a scatter that is 3.7 times worse than with SFD, demon- 
strating that millimeter emission from dust is neither simply related to the Hi 
column density, nor is the dust at one temperature everywhere on the sky. 

Our regression line in Figure § differs significantly from unity, indicating 
that a emissivity is incorrect for the dust. At lower frequencies, the slope 
departs even more strongly from unity. This is an indication that a v'^ emissivity 
is incorrect for the dust, as can be seen from the mean spectrum of large regions 
of the sky (Figure |3|: since the DMR is a differential instrument, we have plotted 
the difference between bright and faint regions of the sky.) 

The SFD prediction can be made using other emissivity profiles by modi- 
fying the exponent, a, in equation ^. An a = 1.5 emissivity profile results in a 
better fit at low frequencies, but ruins the fit at high frequencies. In fact, there 
is no power-law emissivity model that adequately fits the FIRAS data. 
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Figure 2. Comparison of the FIRAS data at 500 GHz with three 
predictors. 



7 




0.0001 L ^^^^^ 

100 1000 
Frequency (GHz) 



Figure 3. Difference spectra from the DMR {u < 100 GHz) and 

DMR {ly > 200 GHz) instruments. These spectra divide the sky 
into three regions based upon temperature, and are renormahzed to 
Aoo/xm = 1 MJy/sr. 
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4. Mult i- Component Dust Model 



Motivated by observations of many different grain species in the ISM, we have 
formulated a general description of a multi-component ISM suitable for mi- 
crowave predictions (Finkbeiner et al. 1999). We capitalize upon the uncertain 
degeneracies between components to construct a multi-component model with 
very few degrees of freedom. As we discussed, any one component of dust might 
be expected to exhibit power-law emissivity in the range 100 <^ y <^ 1000 GHz. 
The complicated behavior of their optical opacities is unimportant. The only 
relevant optical parameter is k*, the effective cross-section to the interstellar 
radiation field (ISRF): 

_ Kk{u)E{v)dv 
Jo E{u)dv 

where E{v) represents the angle- averaged intensity in the ISRF. 

We construct a multi-component model of the ISM, analogous to equation 
^, by summing the emission of these components: 

/ = Y.k fkek{i')Bu{Tk) ^^^^ 

J2k fkekii^o)BuoiTk)Kioo{ak,Tk) 

where fk is the mass fraction of the fc-th component, Tk is the temperature, -f^ioo 
is the DIRBE color-correction factor and /p,ioo is the SFD98 100;um flux. The 
mass fractions are forced to sum to unity, i.e. J2k fk = 1- Each emissivity is a 
power-law with index Ofc. The observationally relevant parameter is the ratio of 
far-IR emissivity to the effective optical opacity, q = e{vQ)/K* . 

Each dust component is therefore described by three global parameters (/fc, 
Qk, OLk) and one parameter that varies with position on the sky, e.g., Ti{x). The 
temperatures of all dust species are coupled by demanding that they all respond 
to the same ISRF (for more details, see Finkbeiner et al. 1999). Whether or 
not the actual components of the dust physically correspond to these compo- 
nents, equation |5| can be thought of as a phenomenological "expansion set" for 
describing the composite dust spectrum. 

Remarkably, we find that only two components of dust are required to 
describe the spectrum everywhere on the sky as seen by FIRAS. As a first guess, 
we tested a two-component model designed to replicate the spectrum in Pollack 
et al. 1994. The prescription in their paper for their best-fit broken power law 
corresponds in our model to ai = 1.5, 02 = 2.6, /i = 0.25, and qi/q2 = 0.61. 
This results in a considerably better fit without fitting any new free parameters. 
Between 800 and 1800 GHz, where the FIRAS signal is very good, the model 
matches the data to approximately 1% everywhere. At lower frequencies, the 
largest deviation is 10%. 

Allowing /i and qi/q2 to float provides even better fits. Our best fit is 
achieved with two components with power-law emissivities ai = 1.68, 02 = 2.78, 
where the former component dominates the emission u ^ 500 GHz. We have 
fit 71% of the sky as observed by FIRAS (4378 pixels) at 123 frequencies, for 
a total of ~ 540, 000 data value comparisons. We also effectively fit for an 
unconstrained zero-point value at each frequency to avoid having the solution 
influenced by uncertainties in the cosmic infrared background. The resulting 
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reduced xt 2.3/DOF, as compared to ~ 30 (!) for the best single-component 
model. Results for a suite of models, including our best fit 2-component model, 
are given in Finkbeiner et al. (1999) along with tests of their robustness. 

5. Discussion 

The thermal emission from Galactic dust can be very successfully predicted 
at millimeter/microwave frequencies using a two-component composition model 
with temperature varying on the sky. We interpret the best fit as a amorphous 
silicate-like component (i/^''^ emissivity, (T) ~ 9.5 K) and a carbonate-like com- 
ponent (z^^'^ emissivity, (T) ~ 16 K). The FIRAS data are also reasonably well 
fit using two v"^ emissivity components, if one of these components is very cold 
((T) 5 K) as proposed by Reach et al. (1995). We find no evidence for a warm 
component {v^ emissivity, (T) ^ 29 K) as proposed by Lagache et al. 1999. 

A full range of Galactic dust models are discussed in Finkbeiner et al. 1999. 
Numerical predictions for thermal emission will be made available at our web 
page, 

tittp: //astro. berkeley.edu/dust , 

or at the FORECAST (FOREgrounds and CMB Anisotropy Scan simulation 
Tools) home page, 

|http: //cf pa. berkeley. edu/f precast . 

By using the IRAS maps as reprocessed by SFD98, these predictions will be 
provided at a resolution of 6fl. 

It is critical that CMBR experimentalists compare their observations with 
valid models for the Galactic dust emission. A "template approach" is often 
carelessly used to compare observations with alleged contaminants, with the 
correlation amplitude indicating the level of contamination. For example, ra- 
dio maps {e.g., the Haslam et al. map at 408 MHz) are used as templates for 
synchrotron emission at much higher frequencies, in spite of the well known 
spectral index variation. 21 cm maps {e.g., Leiden-Dwingeloo) are often used 
as templates for microwave dust emission. This can be a terrible template, as it 
ignores variations in the dust / gas ratio and the well-measured variations in dust 
temperature. This point is demonstrated in Figure where we have correlated 
the Leitch et al. (1997) 14.5 GHz data with possible templates for the dust emis- 
sion. The microwave data shows very little correlation with a 21 cm template, 
or with the even more simple-minded (but still occasionally used) esc |6| model 
for Galactic dust. One might incorrectly conclude from this that there is little 
dust contamination in this CMBR data. In fact, when we correlate the Leitch 
et al. data with a proper prediction at 14.5 GHz, there is evidence for signifi- 
cant contamination. A naive comparison with 100/im IRAS flux also indicates 
a strong correlation with dust - but with the wrong amplitude. We find that 
the dust on the 5 = 88° ring used by Leitch et al. is cold enough to lower the 
100/im emission per dust column by a factor of 2.75 relative to the average high- 
latitude ratio. Therefore, the amount of 14.5 GHz emission per dust column is 
overstated by a factor of ~ 3 in comparisons such as the famous Figure 11 in 
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Figure 4. Correlation of different templates with the Leitch et al. 
1997 CMBR data. A esc |6| model or 21-cm map (from Leiden- 
Dwingeloo) are poor predictors for microwave dust emission. 
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Draine &; Lazarian (1998). Wc hope that our dust spectrum model wih allow 
both a realistic prediction of the thermal (vibrational) dust emission, and also 
facilitate a meaningful comparison between the dust and the mysterious excess 
15-90 GHz emission component that correlates so well with it. 
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